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Abstract 

We calculate elastic electron-deuteron scattering in a chiral effective field theory 
approach for few-nucleon systems based on a modified Weinberg power counting. 
We construct the current operators and the deuteron wave function at next-to- 
leading (NLO) and next-to-next-to-leading (NNLO) order simultaneously within a 
projection formalism. The leading order comprises the impulse approximation of 
photons coupling to point-like nucleons with an anomalous magnetic moment. At 
NLO, we include renormalizations of the single nucleon operators. To this order, no 
unknown parameters enter. At NNLO, one four-nucleon-photon operator appears. 
Its strength can be determined from the deuteron magnetic moment. We obtain not 
only a satisfactory description of the deuteron structure functions and form factors 
measured in electron-deuteron scattering but also find a good convergence for these 
observables. 
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1. A new era of nuclear physics calculations was started by Weinberg [1,2] ap- 
plying effective field theory (EFT) methods and chiral Lagrangians to systems 
of two and more nucleons. Although there are still some (minor) conceptual 
problems in the precise formulation of nuclear effective field theory, it is by 
now established that at very low energies, one can perform very accurate cal- 
culations using a theory of non-relativistic nucleons, whose interactions are 
given in terms of A- nucleon terms (A= 4, 6, . . .) with the pions integrated out, 
the so-called pionless theory. Going to higher energies, the inclusion of pions 
becomes of prime importance and it has been shown that Weinberg's original 
proposal of constructing an irreducible N-nucleon (N= 2,3,.. .) potential and 
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iterating it in a Schrodinger (Lippmann-Schwinger) equation can give a precise 
description of nucleon-nucleon (NN) scattering as well as static and dynamic 
properties of three- and four-nucleon systems, see e.g. [3-5]. The more elegant 
formulation of Kaplan, Savage and Wise (KSW) [6], which allows for power 
counting on the level of the scattering amplitudes, suffers in its present for- 
mulation from the incorrect description of the tensor force, as reflected in the 
non-convergence of the triplet partial waves in np scattering, see [7,8]. For 
a much more detailed discussion of the concepts and applications of nuclear 
EFT, we refer to the recent reviews [9,10]. Of course, many results found in nu- 
clear effective field theory have previously been obtained in more conventional 
meson-exchange approaches [11,12]. These, however, cannot be formulated in 
a truly systematic fashion and cannot be linked simply to the symmetries of 
QCD, as it is the case of the chiral effective field theory employed here | ^'^ | . 

In this letter, we will consider elastic electron-deuteron scattering based on 
a Hamiltonian approach to Weinberg's formulation as developed in [14]. The 
central object will be the (unpolarized) scattering cross section, which in case 
of the deuteron is given in terms of two structure functions. 



with = —Q^ < the invariant momentum transfer squared and 6 is the 
scattering angle in the centre-of-mass frame. Furthermore, we have separated 
the QED (Mott) cross section. Alternatively, one can parameterize the re- 
sponse of the deuteron to an external vector current in terms of three form 
factors, Fc (charge), Fm (magnetic) and Fq (quadrupole). These latter three 
have been the subject of a detailed study in [15], which employed external 
wave functions sewed to the leading order interaction kernel p^ . Our approach 
is similar, only that we construct the current operators and the wave functions 
simultaneously from the same Hamiltonian. Furthermore, we wish to study the 
interplay of nuclear and nucleon dynamics beyond the accuracy than it was 
done in [15]. In addition, we will be able to compare our results directly with 
the pioneering perturbative calculation of [17], which certainly can be consid- 
ered a cornerstone in the development of nuclear effective field theory. Finally, 
this study is to be understood as the first step in a systematic investigation of 
the electromagnetic properties of light nuclei. Since much is known about the 

*'^An illustrative example are isoscalar two-meson-exchange currents. While in EFT 
these have to be small simply due to power counting arguments, in conventional 
meson-exchange models such type of suppression can only be found after often 
tedious calculations, see e.g. [13]. 

^^Earlier related work was performed by Rho and collaborators [16], who worked 
out the current operators to next-to-next-to-leading order (NNLO) and studied in 
particular deuteron disintegration. 
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deuteron from the experimental and the theoretical side, such an investigation 
is clearly needed. More precisely, the deuteron is a very special nucleus due to 
its small binding energy, it thus can be well described by iterated one-pion- 
exchange and some short range physics, as it is done in conventional nuclear 
physics since a long time and has been perfected to high precision. The ap- 
proach presented here is similar but also differs because corrections can be 
calculated in a controlled and systematic fashion. This becomes more impor- 
tant if one goes to more densely bound systems. For a more detailed discussion 
about the usefulness of EFT in the description of deuteron dynamics and the 
relation to conventional approaches, we refer to [15]. 

2. First, we briefly explain the central ideas underlying our calculations. One 

starts from an effective chiral Lagrangian of pious and nucleons, including in 
particular local four-nucleon interactions which describe the short range part 
of the nuclear force, symbolically 

C^eS — C^-K-K + + , (2) 

where each of the terms admits an expansion in small momenta and quark 
(meson) masses. To a given order, one has to include all terms consistent 
with chiral symmetry, parity, charge conjugation and so on. From the effective 
Lagrangian, one derives the two-nucleon potential. This is based on a modified 
Weinberg counting, which is applied to the two-nucleon potential to a certain 
order in small momenta and pion masses. 



with p, p' the nucleon centre-of-mass momenta and the superscript i gives the 
(non-negative) chiral dimension. The power counting underlying this potential 
is based on the considerations presented in [14]. To leading order (LO), this po- 
tential is the sum of one-pion exchange (OPE) (with point-like coupling) and 
of two four-nucleon contact interactions without derivatives. The low-energy 
constants (LECs) accompanying these terms have to be determined by a fit 
to some data, like e.g. the two S-wavc phase shifts in the low-energy region 
(for np). At next-to-leading order (NLO), one has corrections to the OPE, 
the leading order two-pion exchange graphs and seven dimension two four- 
nucleon terms with unknown LECs (for the np system). Finally, at NNLO, one 
has further corrections in the one- and two-pion exchange graphs including 
dimension two pion-nucleon operators. The corresponding LECs can be deter- 
mined from the chiral perturbation theory (CHPT) analysis of pion-nucleon 
scattering (for details, see [4]). The existence of shallow nuclear bound states 
(and large scattering lengths) forces one to perform an additional nonpertur- 
bative resummation. This is done here by obtaining the bound and scattering 




(3) 
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states from the solution of a regularized Lippmann-Schwinger equation. The 
potential has to be understood as regularized, as dictated by the EFT ap- 
proach employed here, i.e. V{p,p') — >• fR{p)V{p,p')f]i{p'), where is a 
regulator function chosen in harmony with the underlying symmetries. Within 
a certain range of cut-off values, the physics should be independent of its pre- 
cise form and value. This range increases as one goes to higher orders, as 
demonstrated explicitly for the np case in [4]. 

3. A powerful method to simultaneously construct wave functions and current 
operators was invented long time ago by Okubo and others [19,20]. The ex- 
tension to meson-exchange currents in few-nucleon systems was pioneered by 
Gari and Hyuga [21,22]. As mentioned before, in [14] it was shown how this 
method could be extended to chiral effective Hamiltonians based on Wein- 
berg's power counting. The corresponding electromagnetic current, denoted 
by throughout, can now be obtained in two ways. One can either gauge 
the effective Hamiltonian or calculate directly the effective current from the 
relation 

i^l^Mi^i) = {(Pjljfm , (4) 

with 10) the projected low-energy state, {i/j) = (1 + A){r]{l + A^A)~^/'^ri)\(l)). 
Here, rj is the projector onto the Fock space with two nucleons and no pions 
and the operator A can be obtained recursively from demanding that the wave 
function orthogonal to r}\ip) decouples. The effective current takes the form 

jf = r){l + A^Ar'/'\l + At) J, (1 + A){1 + A^A)-'/' r] . (5) 

We mention in passing that in this approach one has to include so-called 
wave function reort honor malizat ion diagrams, see e.g. [14] for the two-nucleon 
potential. Their contribution to the effective current is completely cancelled by 
the two-body recoil current (in the non-relativistic limit) [21]. Before giving 
the explicit form of the effective current, it is important to formulate the 
power counting in the presence of external electroweak fields, first worked out 
by Rho [18]. Denoting by Q a small external momentum or a pion mass, any 
QCD matrix-element (in the presence of external fields) takes the form 

A^ = Q^/(^,^) (6) 

where is some regularization scale, g denotes a collection of coupling (low- 
energy) constants and / is a function of order one. The counting index is given 



4 



by 



1 - 2C + 2L + ^ 



di + — + ei 



-2 



(7) 



with C the number of connected pieces in the diagram under consideration, L 
the number of (pion) loops, and P, is the vertex dimension. The latter depends 
on the number of dcrivativcs/pion mass insertions (dj), the number of nucleon 
fields (rii) and the number of external fields (cj), in our case insertions of the 
electromagnetic current. Chiral symmetry demands that the counting index 
is bounded from below, as verified by eq.(7). The terms with C = 2 are called 
one-body terms in the nuclear physics language. They subsume all interactions 
of the photon with either the neutron or the proton (this is also called the 
impulse approximation). For C = 1, obviously both nucleons are involved in 
the interaction, thus one talks of two-body terms (in Weinberg's language, 
these terms are denoted as three-body interactions). From eq.(7), one quickly 
estabhshes that the lowest order terms are of the one-body type and have 
u — —3 (electric photon coupling) and u — —2 (magnetic coupling). These 
comprise the coupling of a photon to point-like nucleons with an anomalous 
magnetic moment and define the leading order. At NLO (and NNLO), we have 
further one-body terms, which are the usual leading one-loop (third order) 
chiral perturbation theory corrections to the nucleon form factors, see [23-25] . 
In addition, we have the leading two-body operators, the celebrated meson- 
exchange (seagull and pionic) currents. Using the leading order vertices as 
given below, these do not contribute to elastic e-d scattering (as it is well 
known). Also, based on this counting, one has no local four-nucleon-photon 
interactions at NLO. The first correction, which is not of one-body nature, 
does appear at NNLO, it is the magnetic four-nucleon-photon interaction 
of [17]. At that order, one has also one-body corrections proportional to the 
magnetic radii of the proton and the neutron. We should point out one subtlety 
with the power counting here. In terms of the wave functions, the current 
operators at NLO and NNLO are only sensitive to the NLO wave functions, 
which can be simply understood from the fact that gauging the corresponding 
operators in the potential does lead to even higher orders in the current. 
Therefore, the results for the electric and the quadrupole form factor at NLO 
and NNLO will coincide, and we denote them as (N)NLO. Formally, one could 
therefore subsume what we call NLO and NNLO here simply as NLO, similar 
to the two orders making up the LO (impulse) contribution. Since this has 
not be done in the existing literature, see e.g. ref. [16], we refrain from doing 
so too. 



We now discuss the effective interaction Hamiltonians including the photon 
field, expressed in terms of the gauge potential and the field strength 
tensor F^,^. Prom these, the effective current follows. Gauging the free and 
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the interaction Hamiltonians, Tio and Tiint, respectively, gives (we only show 
the terms that involve the photon field, all other terms to the order we are 
working are given e.g. in [14]), 



\ rriN imjv / 



(8) 



Kn"^ = ^ ai{-ie^A') t:^ N - eL^ {N^ a, e'^' F,, N) (N^N) , (9) 

with ejv and e^^ the nucleon and pion charge matrix, respectively, and tbat is 
the nucleon mass. We omit terms of higher order in the electromagnetic cou- 
pling since they are not relevant here. Throughout, we work in the Coulomb 
gauge V ■ A = 0. For the gauged free nucleon Hamiltonian, we have only 
given the leading order (dimension one) electric coupling and the (dimension 
two) magnetic one, with Ks the isoscalar anomalous magnetic moment of the 
nucleon, {Kp + Kn)/^ — Ks- These two interactions lead to the terms with 
u = —3 and —2 as discussed before. The dimension three terms contributing 
to the one-body currents at z/ = — 1 can be found in [24]. The pion-photon- 
nucleon interaction in eq.(9) is nothing but the celebrated KroU-Ruderman 
term, with the axial vector coupling measured in neutron /?-decay and 
the weak pion decay constant. As noted before, employing these Hamiltoni- 
ans, the leading exchange currents are proportional to fi x 72 (where 7^ is the 
isospin operator of nucleon i) and thus vanish for isoscalar transitions as it 
is the case for e-d scattering. The leading isoscalar exchange current appears 
only at N^LO, because the only possible NNLO term would consist of a di- 
mension two correction to the KroU-Ruderman vertex. This operator has the 
structure aie^v^q^ (r" + 5^^) (see app. A of [26]) and thus generates only an 
isovector operator. At N^LO, we have contributions from two-pion-exchange 
diagrams with one ■ymrNN-yertex (football and triangle graphs) together with 
dimension three corrections to the pion-nucleon-photon vertex. The last term 
in eq.(9), first considered in [17], is a magnetic photon four-nucleon contact 
interaction, its strength can not be determined from np scattering but can be 
obtained from a fit to the deuteron magnetic moment. We note that there is 
also a NNLO single nucleon correction due to the magnetic radius. We stress 
that such power counting arguments have already been given in [15]. 



4. We now turn to the response of the deuteron to an external electromagnetic 
field. Consider a deuteron state with four-momentum and polarization 

vector e^^, subject to the condition p^e^ = 0. In the deuteron rest frame, 
one conventionally selects = 5f and the corresponding deuteron states are 
denoted by \p, i) . In terms of these and to leading order in the non-relativistic 
expansion, the matrix element of the electromagnetic current is given in terms 
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of the charge, magnetic and quadrupole form factors, 



(p',^|jo|p,j)=e 



if, ^ I j'fe I P, j) = 7, \Fc{Q'^) Sij{p + p')k + FM(Q^){5jkqi - SikQj) 



2md 



(10) 



with q = p' — p, Q = \q\ and is the deuteron mass. These form factors are 
dimensionless and normahzed to the charge, the magnetic moment, fid, and the 
quadrupole moment, Qa, of the deuteron, i.e. -Fc(O) = 1, {e/2m(i)FM{0) = 
and (l/m^)FQ(0) = Qd, with Hd = 0.85741 (e/2miv) and Qd = 0.2859fm2. 
The structure functions A{Q'^) and B{Q'^) defined in eq.(l) are given in terms 
of these form factors via 



B^-ri{l + ri)F, 



2 

M 



Here, r] = (5^/4m^. To disentangle these three form factors, one measurement 
involving polarization is necessary. The analyzing power ^20 has become the 
observable of choice to do that. Since it depends on and the scattering 
angle 9, one also uses the quantity ^20, 



\vFcmFQm + Iv'F'iQ' 



(-20 



V2 



Fhm 



+ 



WFim_ 



(12) 



which only depends on the momentum transfer squared and is independent of 
the magnetic form factor. Note, however, that ^20 is not directly measurable 
(for a detailed discussion, we refer to the recent review [27]). 



5. We first consider the structure functions A and B which can be obtained 
directly from elastic e-d scattering. Since from the study of the nucleon form 
factors in heavy baryon chiral perturbation theory it is known that the third 
order result for the isoscalar form factor starts to deviate from the data at ~ 
0.2 GeV^, we will hmit our range of Q from to 400. . .500 MeV. We work out 
the matrix elements from eqs.(lO) and construct the structure functions by use 
of eqs.(ll) without any approximation. This is different from [17], where the 
kinematical factors like e.g. 7] where also included in the power counting and 
thus in their case there can be no contribution from Fq{Q'^) to A^Q^) at NLO. 
Throughout, we work with an exponential regulator fuip) = exp(— p'^/A*^) 
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with a cut-off A = 600 MeV (similar to what was done in the simultaneous fit 
to np and pp scattering phases, see [28]). None of the results presented here 
depends notably on the cut-off within its allowed bounds. For completeness, 
we give the corresponding LECs. At LO, we have C^si = 0.0427 ■ lO^GeV"^ 
and at NLO, we get C^si = -0.0368-10^ GeV'^, C^si = 0.186-10^ GeV"^ and 
Cssi-sDi = —0.190 ■ lO^GeV"^. These numbers are in good agreement with 
what has been found before using such type of regulator, see e.g. ref. [29]. 




0.1 0.2 0.3 0.4 



Q [GeV/c] 

Fig. 1. The structure function A{Q'^). The dashed and the solid line give the LO, 
the (N)NLO result, respectively, and the data are from [30], [31]. As explained in 
the text, the NLO and NNLO results for A{Q^) can not be distinguished on the 
scale of the figure. 

Let us discuss A^Q"^) as given in fig. 1. Note that we plot the structure function 
versus Q = \fQ^ (not versus Q^) to facilitate the comparison with ref. [17]. 
The leading order curve already gives a fair description of the data, the NLO 
and NNLO result is visibly improved but still above the data. Note that the 
difference between NLO and NNLO stems from a Foldy-type contribution 
proportional to the magnetic radius, which is very small, and an induced con- 
tribution from the magnetic photon-four-nucleon interaction. Both of these 
are small and further suppressed by a factor of r], so that the difference between 
the two curves can not be seen on the scale of fig. 1. Since there is no exchange 
current contribution, the improvement going from LO to NLO is due to the 
single nucleon radius terms which appear at this order. This can be seen by 
considering the LO wave function together with the NLO current operators. 
At yet higher orders, however, there should be some further improvement due 
to the wave functions as shown in [15]. One should consider such a sensitiv- 
ity to single nucleon observables not as an unwanted complication but rather 
conclude that indeed isoscalar quantities (alas neutron properties) can be in- 
ferred from nuclear targets with some precision. We consider this interplay of 
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chiral nucleoli dynamics and nuclear EFT as one of the important ingredients 
in the whole approach. Note further that our (N)NLO result is similar to the 
one of [17]. It should also be noted that for the range of momentum transfer 
considered here, the NLO corrections are sizeably smaller than in the KSW 
calculation, pointing towards a better convergence. 

Next, we turn our attention to B{Q'^), which is essentially the response to the 
magnetic photon, see fig. 2. Again, the LO result is visibly better than the one 
obtained in [17], whereas the NLO curves are comparable. However, in contrast 
to the work of [17], our prediction for B{Q'^) is parameter-free, whereas the 
one in the KSW scheme has one free parameter from the four-nucleon-photon 
interaction (due to the different counting). Also, the NLO correction is fairly 
small for the momentum range considered here, indicating convergence. At 
NNLO, we have to pin down the LEG L2 of the magnetic photon four-nucleon 
interaction, which is done by fitting the magnetic moment of the deuteron. 
The resulting curve for B{Q^) is further improved, see fig. 2. 
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0.1 0.2 0.3 0.4 

Q [GeV/c] 

Fig. 2. The structure function B(Q^). The dashed, the dot-dashed and the solid Une 
give the LO, the NLO and the NNLO result, respectively. The new data are from 
the compilation of Sick [32] (filled diamonds) and the older data (open circles) are 
from [33], [34], [35]. 

Next, we discuss the charge, magnetic and quadrupole form factors. The cor- 
responding normalizations (static properties) are collected in table 1 for LO 
and NLO. The only difference at NNLO is the magnetic moment, which is 
exactly reproduced (since it is used as input to fix L2). The corresponding 
value for the the LEG L2 is L2 = 2.68 ■ lO^^GeV^^, which is rather small. 
This can be traced back to the fact that the NLO result for is already 
within 1% of the empirical value. Note that in contrast to what was done in 
ref.[4], we have fine tuned the LEGs in the deuteron channel to reproduce the 
binding energy to four digits. The only difference at NNLO to NLO is the 



9 



magnetic moment, therefore the NNLO numbers are not given. Notice that 
with the fine tuned binding energy, the prediction for the quadrupole moment 
is visibly improved as compared to ref. [4] and closer to the data than in all 
modern high-precision potentials. 

Table 1 

Static properties at LO and NLO. Here, Ed, and Qd denote the deuteron binding 
energy, its magnetic and its quadrupole moment. 





LO 


NLO 


Exp. 


\Ed\ [MeV] 


2.224 


2.224 


2.22456612(12) 




0.828 


0.852 


0.8574382284(94) 


Qd [fm2] 


0.265 


0.276 


0.2859(3) 



In all cases, the NLO/NNLO corrections are small and one finds a decent 
description of the three form factors, see fig. 3. The largest discrepancies at 
higher momentum transfers are found in the charge form factor, as refiected in 
the results for A(Q^) shown above. Similar results have been obtained in [15], 
however, a direct comparison with data was not given in that paper. 




Fig. 3. Electric, Fc{Q'^), magnetic, Fm{Q^), and quadrupole, Fq{Q'^), form factors 
of the deuteron, in order. The dashed, the dot-dashed and the soUd Unes give the 
LO, the NLO (only shown for Fm) and the NNLO results, in order, and the data 
are from the analysis of Sick [32]. As explained in the text, NLO and NNLO results 
coincide for the electric and the quadrupole form factor. 



From the electric and the quadrupole form factor, we can construct t20) as 
shown in fig. 4. Again, the deviations of our LO and (N)NLO predictions for 
Fc(Q^) from the data are responsible for the too small magnitude of t2o at 
Q ~ 0.4 GeV. 
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Fig. 4. Analyzing power t20 at LO (dashed) and (N)NLO (solid line). The data are 
from the recent compilation [36] . 

6. We have analyzed electro n-deuteron scattering in the framework of a chiral 
effective field theory for few-nucleon systems at next-to-leading and next-to- 
next-to-leading order. At NLO, no meson-exchange currents or four-nucleon- 
photon operators contribute. At NNLO, only one magnetic photon-four-nucleon 
operator appears, the corresponding coupling constant can be determined from 
a fit to the deuteron magnetic moment. As stressed, the NLO and NNLO pre- 
dictions for the electric and the quadrupole form factors coincide at NLO and 
NNLO because one is not yet sensitive to the NNLO wave function corrections. 
In particular, we have discussed the interplay between the single nucleon dy- 
namics as encoded in the nucleon form factors and the nuclear dynamics. As 
already stressed in [15], the accuracy of the description of the nucleon form 
factors limits the applicability of the effective field theory approach to the 
deuteron structure to momentum transfer of about (5 — 0.4 GeV. Thus an 
improved description of these single nucleon observables has to be obtained 
to extend these considerations consistently to higher photon virtualities. For 
a step in this directions, see [25]. We stress again that we consider it impor- 
tant to consistently describe the single as well as the few-nucleon sector. The 
results presented here extend the ones of [15] to the first non-trivial order 
for the single nucleon form factors, i.e. are sensitive to the structure of the 
nucleon. Obviously, as next steps one has to consider N^LO corrections to the 
process discussed here (since only at that order sufficiently many additional 
contributions appear, like e.g. meson-exchange currents) as well as electron 
scattering of three- and four-body systems, using the wave functions of [5]. 
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